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Summary 


This manual describes the use of the computer code ICAN/DAMP (Integrated 
Composite Analyzer with Damping Analysis Capabilities) for the prediction of damping 
in polymer-matrix composites. The code is written in FORTRAN 77 and is a version of 
ICAN (Integrated Composite ANalyzer) computer program. The code incorporates a 
new module for synthesizing the materials damping from micromechanics to laminate 
level. Explicit micromechanics equations based on hysteretic damping are programmed 
relating the on-axis damping capacities to the fiber and matrix properties and fiber 
volume ratio. The damping capacities of unidirectional composites subjected to off-axis 
loading are synthesized from on-axis damping values. The hygrothermal effect on the 
damping performance of unidirectional composites caused by temperature and moisture 
variation is modeled along with the damping contributions from interfacial friction 
between broken fibers and matrix. The temperature rise in continuously vibrating 
composite plies and composite laminates is also estimated. The ICAN/DAMP users 
manual provides descriptions of the damping analysis module’s functions, structure, input 
requirements, output interpretation, and execution requirements. It only addresses the 
changes required to conduct the damping analysis and should be used in conjunction with 
the ‘Second Generation Integrated Composite Analyzer (ICAN) Computer Code” users 
manual. 


Introduction 


CAN/DAMP (Damping Analysis of Composite Laminates), a version of the 
CAN (refs. 1 and 2) computer code, incorporates a new module for predicting damping 
in polymer-matrix composites. The new module allows ICAN to synthesize material 
damping at all scale levels of the composite— from micromechanics to laminate level— 

thus, providing the required data for simulating the damped vibrational characteristics of 
composite structures. 

This manual provides a step-by-step guide to using the new damping module and 
only addresses those changes required to conduct the damping analysis. Therefore, users 
should be familiar with the input requirements and analysis capabilities of the standard 
ICAN code as described in references 1 and 2. The ICAN/DAMP users manual de- 
scribes the damping analysis module’s functions, structure, input requirements, output 

podtuam 011 ’ and execution re< l u * r ements. The ICAN/DAMP code is compatible with 
t OR l RAN 77 standards; however, the system requirements of various hardware environ- 
ments may differ, and the user should consult computer support personnel in the event of 
insta ation problems. Although the results of the damping module have been sufficiently 
ested, the user is solely responsible for the interpretation and application of the obtained 

results. Questions regarding the damping module should be directed to the authors of 
the manual. 


Damping Module Overview 


Background 

A brief overview of the theoretical background of the damping theory for compos- 
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ites and Tiber composite laminates is discussed in this section. Detailed descriptions of 
the damping theory, including equations and their application in the analysis, are 
provided in references 3 and 4 . Although some micromechanical damping and laminate 
damping equations are included in appendix l 1 , we strongly suggest that the user read 
references 3 and 4 for further understanding of the integrated damping theory. 

The damping of composite mateiials is anisotropic and depends on several micro- 
mechanical, laminate, and structural parameters. The damping of on-axis composites 
(unidirectional composites loaded along the material axes as shown in Tig. 1 (a)) is 
expressed in terms of six specific damping capacities (SDC), of which only four are 
independent. These SDC’s are defined as the dissipated strain energy of the material 
over the maximum stored strain energy during one vibration cycle. These SDC s are 
longitudinal (^ iu ), transverse in-plane normal {^ m ), in -P lane sliear (^12)’ through-the- 
ft bers normal (^733)5 through-the-fibers shear (V 1 /^)? an< l out-of-plane shear (^23)' 
on-axis SDC’s are related to the mechanical properties of the constituents (fiber and 
matrix), the fiber volume ratio (FVR), temperature, moisture, and existing rnicrodamage 
(percentage of broken fibers). In addition to the previous parameters, the SDC’s of off- 
axis composites (unidirectional composites loaded at an angle with respect to the fiber 
direction) depend on the fiber orientation angle (see fig. 1 (b)). 

The damping of composite laminates includes the damping contributions of all 
plies. In addition to the previously described parameters affecting 011-axis SDC’s, the 
damping of composite laminates depends on ply angles, ply thicknesses, and the overall 
ply stacking sequence. A unique feature of the encoded theory is its capacity to predict 
the damping of general laminates as it includes the damping effects of extension, flexure, 
and extension-flexure coupling. These different sources of laminate damping are best de- 
scribed by the extensional [Aq], flexural [Dp], and extension-flexure coupling [Cp] 
damping matrices. The overall damping capacity of a composite laminate is related to 
these three damping matrices and the local deformation expressed by midplane strains 
and curvatures. 

Increased temperatures may be observed in vibrating composite laminates as a 
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Figure 1 .—Unidirectional composites. 
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result of strain-energy dissipation, including the friction effects. It is well known that 
damping represents the strain/kinetic energy losses that are dissipated in the composite 
material in the form of heat. For extended periods of vibration, it is possible that a 
significant temperature rise may occur, which affects such mechanical properties of the 
laminate as stiffness and strength. It is desirable, then, to predict both the temperature 
rise in composite laminates during prolonged vibrations, and its effect on the mechanical 
properties. The damping module provides this capability by calculating the temperature 
profiles developed through the thickness of the vibrating laminate. 


Module Functions 

The present module was developed for analyzing the damping capacity of 
unidirectional composites and/or composite laminates. The modeling of damping is 
performed in two levels: (1) the micromechanics level, where the SDC’s of unidirectional 
on-axis composites and/or off-axis composites are synthesized, and (2) the laminate level, 
where the damping capacity of composite laminates is calculated. 

The input information for the damping analysis requires additional card groups to 
control the output of the module and to provide additional input parameters. These 
damping analysis card groups are part of the input data file of the main ICAN code. 
Detailed descriptions of the damping analysis card groups are provided in the “Input 
Data section of this manual. In addition to the damping analysis card groups, the 
module uses the standard input parameters previously available through the ICAN input 
file. These parameters are 

(t) Data defining the laminate configuration and environmental conditions such 
as ply angles, ply thicknesses, ply composite materials, ply cure temperature, ply use 
temperature, and ply moisture content 

(2) Definition of the basic composite material systems in terms of fiber volume 
ratios (FVR’s), fiber material, and matrix material 

(3) Mechanical loads applied on the laminate. 

Details of these parameters are addressed in the input data sections of the ICAN user’s 
manuals (refs. 1 and 2). 

Ihe properties of the fibers and matrix for each composite material, including 
specific damping capacities, are provided in a separate data base of material properties 
(see appendix C). The user can easily incorporate additional constituent materials based 
on need. A detailed description of the material data base is included in the “Data Bank” 
section of this manual. 

The calculated SDC s of the constituent (fiber and matrix) materials, the compos- 
ite plies, and the composite laminate, and the resultant temperature profiles through the 
thickness of the laminate are printed in a separate section of the output file. Details 
about the damping analysis output are given in the output section of this manual. Also, 
a brief description of the damping module’s subroutines is given in appendix G. 


Input Data 


Input to the damping module is provided by two files: (l) the ICAN/DAMP 
input data file and (2) the constituent materials data base. Typical examples of these 
files are shown in appendices B and C, respectively. 
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Input Data Requirements 

The following input parameters are used by the damping module. This informa- 
tion should be known by the user before executing the code. These parameters are 

(1) Data defining the laminate configuration and environmental conditions such 
as ply angles, ply thicknesses, ply composite materials, ply cure temperature, ply use 

temperature, and ply moisture content . 

(2) Definition of the basic composite material systems m terms ot fiber volume 

ratios (FVR’s), fiber material, and matrix material 

(3) Mechanical loads applied on the laminate t .... 

(4) Properties of the fibers and matrix for each composite material, including 

specific damping capacities. „ . , . 

Optional input data parameters may be needed for the damping analysis as ex- 
plained in the following subsections. Optional input may include 

(1) Amplitudes of cyclic mechanical loads (or laminate strains) applied on the 
laminate and the corresponding forcing frequency. The mechanical loads will affect only 
the equivalent laminate damping, temperature rise, and friction results. 

(2) Data describing the percentage of broken fibers and friction coefficients il 
friction damping contributions due to broken fibers are requested. 


ICAN/DAMP Input File 

The data cards of the input file are summarized in table I. The ICAN/DAMP 
input file is divided in two data sections. The first section is the standard ICAN input 
described in reference 2 and table I. Again, it is suggested that the user consult ref- 
erence 2 for descriptions of other cards included in the input data file. The second data 
section incorporates the damping analysis card groups as shown in appendix D. J he 
following is a detailed description of the damping analysis card groups with examples 

taken from appendix B. , , . .. 

Damping option card group V. - This set of three cards controls various options 

of the damping module: 


Damping option value 

Optional description 

-1- -2- 

-3- 

-4- 

-5- -6- -7- -B- -9- 

1 

frtctn 

— 1 

0 IF: 

1 - 

— - 1 1 1 1 
Include Friction Effects 

TEMPRTS 

0 

1 - 

Temp ' s from Individual Ply Heat Generation 

DAMPOUT 

0 

1 - 

Suspends Damping's Output Printout 


Each card ha, a forma!. (8X.I8) beginning with “FRICTN", “TEMPMS” and “DAMP- 
OUT’ as their first entry, followed by the damping option value under field -1-. lie 
default values of the integer flags are 0. These nag values have the following functions: 
in FRICTN = 1: The friction effects on damping due to broken fibers is mclud 
ed; otherwise, it is neglected (FRICTN - 0). (The default value of FRICTN is 0) 

(2) TEMPRIS = 1: The temperature rise is obtained from calculations ot tne 
dissipated strain energy in each ply; otherwise, the temperature profiles are calculated 
from the average laminate dissipated energy default value of TLMl lvlo - U. 

(3) DAMPOUT = 1: The printout of the damping output is suspended. 
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TABLE I. — SUMMARY OF DETAILS FOR PREPARING ICAN/DAMP INPUT DATA CARDS 


Card 

group 

Identifica- 

tion 

Code symbol 

Number 

of 

entries 

List of en- 
tries 

Card 

field 

col- 

umns 

Format 

Comments 

1 

Title card 

TITLE 

80 

Alphabetic 

character 

1 to 80 

(10A8) 


2 

Booleans 

COMSAT 

1 


1 to 
16 

(A8,L8) 

T(true) if laminate 
analysis is desired; 
otherwise F (false) 



CSANB 

1 


ft 


T (true) if symmetry 
exists; otherwise F 
(false) 



BIDE 

1 


ft 


T (true) if inlerply 
contributions are de- 
sired; otherwise F 
(false) 



RINDV 

1 


ft 


T (true) if disps. are 
input; otherwise F 
(false) 



NONUDF 

1 


ft 


T (true) if Poisson’s 
ratio differences 
chart is not desired; 
otherwise F (false) 



DEFECT 

1 


ft 


T (true) if durability 
analysis with defect 
is desired; otherwise 
F (false) 

3 

Ply (ply 
desired) 

JNPljIPl, TU, TCU, 
DELM, TETA, 
TI1CKNS 

7 

b ji T u , T cu , 
AM, 0j, tj 

1 to 64 

(A8,2I8, 

5F8.3) 

Ply layup, temper- 
ature and moisture 
conditions 

4 

MATCRD 

(material) 

IPl, CODES (1,1, J), 
CODES (l,2,J), VFP, 
VVP, CODES (2,1,J), 
CODES (2,2, J), VFS, 
VVS, VSC 

8 

Primary 
composite 
code naine9 
for fiber and 
matrix, k f , 
k v ; secondary 
composite 
code names 
for fiber and 
matrix k 9C , 

kf' k v 

1 to 72 

(A8,I8, 

2A4, 2E8.2, 
2A4, 3E8.2) 

Descriptions of the 
material system to 
be used 
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TABLE I.— Concluded. 



Card 

group 

Identification 

Code symbol 

Number 

of 

entries 

List of entries 

Card 

field 

columns 

Format 

Comments 


5 

FRICTN 

IOPT(l) 

1 

0 (Default 
yalue) or 1 

1 to 10 

(8X,I8) 

Integer flag to 
include friction 
due to broken 
fibers 

I 

i 


TEMPRIS i 

IOPT(2) 

1 

w 

0 

0 

Integer flag to 
include tempera- 
ture rise from 
individual ply 
heat generation 

i 

■ 


DAMPOUT 

IOPT(3) 

1 

K 

0 

0 

Integer flag to 
suspend the 
damping’s output 
printout 


6 

FREQ 

FREQ 

1 

f 

1 to 16 

(8X.F8.3) 

Forcing frequency 
in herts 

= 

7 

STRAINS 

E0(1),E0(2), 

E0(3) 

3 

* x »*y>7xy 

1 to 32 

(8X.3F8.3) 

User-defined mid- 
plane strains for 
the damping anal- 
ysis 



CURVTUR 

CURV(l), 

CURV(2), 

CURV(3) 

3 


1 to 32 

0 

User-defined cur- 
vatures for the 
damping analysis 


S 

PLY (Same 
mnemonic as 
in card group 

3) 

IL,MU, KFB 

3 

Nn, 

1 to 32 

(8X, 18, 
2F8.3) 

Ply identification 
number, fiber/ 
matrix friction 
coefficient and 
fraction of broken 
fibers in a ply 
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PMEMB 

NX, NY, NXY, 
THCS 

4 

N x >N y iN xy ,O cs 

1 to 40 

(A8, 4E8.4) 

Membrane loads 
and angle to stru- 
ctural axis 



PBEND 

MX, MY, MXY 

3 

Mx.My.Mxy 

1 to 32 

0 

Bending loads 



PTRAN 

DMX, DMY, 
PU,PL 

4 

P u' P l 

1 to 40 

0 

Traverse loads 
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Cyclic loads 
cards 

CNXX, 

CNYY, CNXY 

CMXX, 

CMYY, 

CMXY 

4 

Upper limit, 
lower limit, 
No. eye. and 

Pi 

1 to 40 

(A8.4E8.4) 

Description of 
cyclic loads 


11 

PRINT (out- 
put option 
cards) 

IDECIIO, 

INPTSUM, 

1 


1 to 16 

(2A8) 

Select output 
items 




ALL 

- 






i 
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Forcing frequency: card group VI. — The card in this group has a format 
(8X,F8.3) beginning witii the mnemonic “FREQ” as the first entry and is followed by 
the value of the forcing frequency f (in hertz) as its second entry: 


Frequency value 

- 1 - - 2 - 

1 1 

FREQ 50.0 


Note that the laminate loads in the input are considered in the damping analysis as 
amplitudes of cyclic loads of frequency f, which are used in the calculation of the 
equivalent laminate SDC. Even though this frequency value is needed only for the 
calculation of the temperature profiles, the presence of a forcing frequency card is re- 
quired in the input file. For this example, the forcing frequency is equal to 50.0 Hz. 

Deformation: card group VII. — These two cards provide the option to calculate 
the laminate damping and temperature profiles for a user-defined harmonic laminate 
deformation (expressed by the midplane strain and curvature amplitudes) of frequency f: 





^xy 

-1- 

1 

-2- 

i 

-3- 

i 

1 

1 

STRAINS 

| 

0.0 

- 1 “ 
0.0 

1 

0.0 










- 1 - 
1 

-2- 

| 

-3 - 

| 

-4- 

I 

1 

CURVTUR 

1 

0,0 

1 - 

0.0 

1 

0.0 


In the default option, the laminate deformation is calculated from the sinusoidal load of 
frequency f defined in the loads card group. The default option is assumed when all 
entries in the previous two cards are set equal to 0.0 (as in this example). The format of 
each card is (8X,3F8.3). 

The entries for the midplane strains begin with the mnemonic “STRAINS” as the 
first entry under field -1-, followed by three real entries consisting of the prescribed 
values of e x , s y , and 7 xy , respectively. Each real field is entered in format F8.3. The 
entries for the curvature card begin with the mnemonic “CURVTUR,” followed by the 
values of « x , /c y , and /c xy , entered also in format F8.3. 

Friction data: card group VIII. — The last group of cards in the damping 
analysis input file consists of friction data for each ply used in the calculation of friction 
damping due to broken fibers. 
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Ply id 


"lb 

-1- 

! 

-2- 

1 

-3- 

1 

-4- 

1 

1 

PLY 

1 

1 

1 

0.30 

l 

0.20 

PLY 

2 

0.30 

0.20 

PLY 

3 

0.3C 

0.20 

PLY 

4 

0.30 

0.20 

PLY 

5 

0.30 

0.20 

PLY 

6 

0.30 

0.20 

PLY 

7 

0.30 

0.20 

PLY 

8 

0.30 

0.20 


The format of each card is (8X,I8,2F8.3) and consists of the ply identification number 
(second entry), the friction coefficient between matrix and fibers (p) for the ply (third 
entry), and the fraction of broken to total fibers (N^) in the ply (fourth entry). The 
friction effects are accounted for on the laminate SDC, and the programmed formulation 
is valid for uniaxial cyclic stress only. It is assumed that the modulus of the composite 
does not change with the percentage of broken fibers. Therefore, the stiffness of the 
ply /laminate will not be affected. The ply identification numbers should be in sequential 
order starting from 1, and the values of p and should be specified for each ply. 

The number of plies in this card group should match the number of plies specified in the 
laminate configuration card group. The values for p and used in the example are 

p = 0.30 and Nfj, = 0.20. This option is activated through the damping option card 
group by setting FRICTN = 1. 


Data Bank 


A sample data base of constituent properties is provided with the source code to 
facilitate the user’s initial interaction with the program. After getting acquainted with 
the damping module, the user may want to expand the data bank of constituent proper- 
ties based on his/her specific applications. Among other things, the resident data bank 
reduces the burden of repeatedly preparing input for the fiber and matrix properties. 

The data b ank is structured in the same way as described in the “Second Generation 
ICAN Computer Code Users Manual” (ref. 2). This section is written with the assump- 
tion that the user is already familiar with ICAN’s revised data bank structure and fea- 
tures; therefore, only the additional cards will be discussed. The additions to this data 
bank consist of the fiber and matrix damping capacities cards necessary for the damping 
analysis. These additional cards are described below. 

Fiber material properties. — There are four additional cards that were added to 
the block of fiber properties, corresponding to their specific damping capacities in the 
longitudinal (V’fn)i shear (V T fl2 ), and transverse directions (V’f22» ^23)- lh ese f° ur car ds 
are included after the property fiber compressive strength, and their numerical values are 
entered starting from column 42 in format E9.3. This will increase the total number of 
cards in the fiber properties “block” to 24. An example of the fiber properties with the 
additional fiber damping capacities cards (identified by bold print) is shown in the 
following table: 
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Description 


Symbol 


Value 


Units 


-l- 


- 2 - 


“f- 


-3- 


-I- 


- 6 - 


-7- 


- 8 - 


HIGH MODULUS SURFACE TREATED FIBER. 


HMSF 
$ 

$ 

$ 

Number of fibers per end 
Filament equivalent diameter 
Weight density 
Normal moduli (11) 

Normal moduli (22) 

Poisson's ratio (12) 
Poisson's ratio (23) 

Shear moduli (12) 

Shear moduli (23) 

Thermal expansion coef. 
Thermal expansion coef. 

Heat conductivity (11) 

Heat conductivity (22) 

Heat capacity 
Fiber tensile strength 
Fiber compressive strength 
Fiber damping capacity (11) 
Fiber damping capacity (22) 
Fiber damping capacity (12) 
Fiber damping capacity (23) 


( 11 ) 

( 22 ) 


Nf 

10000 

number 

df 

0 . 300E-03 

Inches 

Rhof 

0. 703E-01 

lb/in* *3 

Efll 

0 . 550E *08 

psi 

Ef 22 

0 . 900E * 06 

psi 

Muf 12 

0. 200E + 00 

non -dim 

Nuf23 

0 . 250E + 00 

non-dim 

Gf 12 

0. 110E+07 

psi 

Gf23 

0. 700E+06 

psi 

Alfafll 

- - 550E-06 

in/in/F 

Alfaf22 

0 . 560E-05 

in/in/F 

Kfll 

0. 403E+01 

BTU-in/hr/in** 2 /F 

Kf22 

0 . 403E+00 

BTU-in/hr/in* * 2 /F 

Cf 

0. 170E+00 

BTU/lb/F 

SfT 

0. 280E+06 

psi 

SfC 

0. 200E+06 

psi 

Zotafll 

0.400E-02 

non-dim 

Zetaf22 

0. 400E-02 

non-dim 

Zetafl2 

0. 40OE-O2 

non-dim 

Zetaf 23 

0. 400E-02 

non-dim 


Matrix material properties. - Two additional cards were included in the block of 
matrix propert.es. These cards correspond to the matrix normal (V- ,) and shear U ) 
specific damping capacities, which are included after the property glass transition 
temperature. Again, the numerical values are entered starting from column 42 in format 
*,9.3. The total number of cards in this block of matrix properties increased to 21 An 
example depicting the matrix properties available in the data bank is shown in the list 
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Description 


Symbol 


Value 


Units 


-1- 


- 2 - 


-5- 


- 8 - 


1 1 1 

0210 DX210 RESIN - 

$ 

$ 

$ 

Weight density 
Normal modulus 
Poisson's ratio 
Thermal expansion coef- 
Matrix heat conductivity 
Heat capacity 
Matrix tensile strength 
Matrix compressive strength 
Matrix shear strength 
Allowable tensile strain 
Allowable compr. strain 
Allowable shear strain 
Allowable torsional strain 
Void heat conductivity 
Glass transition temperature 
Matrix damp, capacity-normal 
Matrix damp. capaclty-Shear 


Rhom 

Em 

Hum 

A1 fa m 

Kill 

Cm 

SmT 

SmC 

SmS 

eps mT 
eps mC 
eps mS 
eps mTOR 
kv 

Tgdr 

Zeta mn 
Z«ta ms 


0. 4 40E-01 

lb/ in* * 3 

0 . 500E + 06 

psi 

0 . 350E * 00 

non-dim 

0 . 360E-04 

in/in/F 

0. 868E-02 

BTU-in/hr/in* A 2/F 

0 . 250E + 00 

BTU/lb/F 

0. 150E+05 

psi 

0. 350E + 05 

psi 

0. 130E+05 

psi 

0 . 200E-01 

in/in 

0 . 500E-01 

in/in 

0 . 35UE-01 

in/in 

0. 350E-01 

in/in 

0. 225E + 00 

BTU-in/hr/in**2/F 

0 . 4 2 0E f 0 3 

F 

0.103E+00 

non-dim 

0 . 118E+00 

non -dim 


Output 


Tlie output is organized in Tour sections describing (1) the : on-ax, s damping SDC s 
,f earl, composite material, (2) the off-axis SDC’s of each ply (3) the eqnr.alent 
lamping characteristics of the laminate, and ( 4 ) the calculated temperature profiles A 
.epical output sample is provided in appendix D. More analytical descriptions of the 
lamping output section are given in the following sections: 


i 

I 


t 


i 


On-axis Damping 


This portion of the output provides details of the 
material under a use temperature and moisture content, 
fibers, matrix, and the composite material are obtained: 


on-axis SDC’s of each composite 
The following SDC’s of the 


Specific Damping Capacities 

Fiber 

Matrix 

Composite 

ZF1 1 = ^ fll 
ZF12 = tp n2 
ZF13 = $fxz 
ZF22 = ^f22 
ZF23 = V^f23 
ZF33 = ^f33 

ZMN = ^ nn 
ZMS = 

ZL1 1 — in 
ZL12 = ii 2 
ZL1.3 = 

ZL2i — V^/21 
ZL22 = if) J22 
ZL23 = ^/23 
ZL31 = 0 J31 
ZL32 = !^/32 
ZL33 = 
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OfT-axis Damping 


lip , T f e off - a ™ dam P«ng matrix and the axial damping with respect to the structural 
axes x,y) of each ply are printed out. The axial SDO’s represent the global damping of 

(lef 3) hCn 13 SUbjeCted t0 un,axial cyclic normal and shear stress in the x-y plane 


Damping Matrix Terms 

Axial Damping 

ZC (M) = tp 
ZC(i,2) = 0cxy 
ZC(1,3) = ^ cxs 
ZC(2,1) = ^ 
ZC(2,2) = A 
ZC(2,3) = * 

ZC (3.«) = ^ csx 

ZC(3,2) = # 
ZC(3,3) = ^ css 

ZCAXLX = 0 

„ *x, axial 

ZCAXLY = 4, 

_ ’y, axial 

ZCAXLS = tl> . , 

axial 


Laminate Damping 

ThlS Se ^| 0n i 0f U,e out P ut Provides data related to the total damping capacity of 
e am mate. I he laminate’s midplane strain and curvature are first printed, followed 
y ie equiva ent SDC of the laminate, the dissipated strain energy per unit area and 
the maximum stored strain energy per unit area. These values ar e aS related lo Uie 
previously mentioned generalized laminate strains. 


Laminate mechanical 
strains 

Midplane 

strains 

Curvature 

E01 = e 
E02 = fy 
E03 = 7 

y 

KI = « 
K2 = * y 
K3 = * xy 


SDC of laminate for above strains 

0 = 6 W/ W 

Dissipated strain energy 

6 W 

Stored strain energy 

W 


the lamWi 'T""*'' ‘'’‘"’’’T M I D d). [C D I representing, respectively, 

e laminate damping capacities m extension, flexure, and extension-flexure action are 
also provided m this output section. 


11 





Extensional Da mping [ ] 

ACD (1,1) • • * ACD (1,3) 

ACD (2,1) ACD (2,3) 

ACD (3,1) • ■ - ACD (3,3) 


Laminate Damping Matrices 


Coupl i ng Damping (Cl 
CCD (1,1) ... CCD (1,3) 
CCD( 2,1) ... CCD( 2,3) 
CCD( 3,1) ... CCD (3, 3) 


Flexural Damping [D J 
DCD (1,1) DCD (1,3) 

DCD ( 2 , 1) ... DCD (2, 3) 

DCD (3,1) ... DCD( 3,3) 


Temperature Profiles 

The developed temperature differentials in each ply of the vibrating laminate due 
to the dissipated strain energy are provided in this output portion e lnitia einpera 
tore the average temperature differential, and the temperature differential at the bottom 
and ’top of the ply are printed for each ply. Currently, tl.e heat dissipation is being done 
by heat conduction through the material. The air convection rate is neglected at this 

lime. 


Execution and General Remarks 

ICAN/DAMP has been compiled using a variety of Fortran 77 compilers on 
hardware ranging from microcomputers to supercomputers with relative ease; therefore, 
it is not foreseen that the user would have any great difficulty compi in 8 ie co 
ICAN/DAMP was primarily compiled and executed on the VM/SP operating sys ,em, 
residing in an Amdahl 5870 mainframe computer at NASA Lewis Research Center. The 
same REXX EXEC files that compile and run the ICAN code under the VM system are 
applicable to the ICAN/DAMP code. The input and output files are associated to the 

following I/O units: 

(1) ICAN/DAMP input file I/O unit 5 

(2) Data bank file I/O unit 8 

(3) ICAN/DAMP output file I/O unit G 

These JCL (job control language) Hies are si, own in appendix E. Also a DCL (digital 
command language) command procedures We is provided to execute tl.e code a 
VAX/VMS environment. 


Demonstration Problems 

In addition to the input file, the following demonstration problems are also 
provided in order to aid the user in becoming familiar with some of the code s features 

and capabilities. 


Problem 1 

This problem illustrates the interfacial frictional damping contribution caused by broken 
fibers Tn an off-axis composite. A two-ply, 0.02-in,lhick continuous-fiber compose is 
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subjected to an off-axis, axial stress amplitude <r cxx of 20 ksi (N = 400 lb/in.). The 
composite is made of HMSF (high modulus surface treated graphite fibers) and IMIIS 
(intermediate modulus high strength) matrix with a 0.50 fiber volume ratio. The 
coefficient of friction between fiber and matrix was assumed to be 0.30; the percentage of 
broken fibers per ply to equal 20 percent; and the fiber orientation angle Oj to be 20 . 

The input data set for this problem is given below: 


$ This problem illustrates the interfacial frictional damping 
$ contribution in the case of off-axis loading. 

$ 

I CAN/DAMP: Demonstration Problem No. 1 

COMSAT T 

CSANB F 

BIDE F 

RINDV F 

NONUDF T 

DEFECT F 

$ The composite is 0.02 inches thick and the fiber 
$ angle is equal to 20 degrees. 

PLY 1 1 70.0 70.0 0.0 

PLY 2 1 70.0 70.0 0.0 

$ HMSF fiber in Intermediate Modulus High Strength 
MATCPD 1HMSFIMHS 0.50 0.00T300IMHS 

$ Friction effect is considered in this problem. 


orientation 


20.0 

20.0 

matrix 

0.00 


0.01 
0.01 
( IMHS ) . 
0.00 


0.00 


FRICTN 

1 

IF: 1 

- Include Friction Effects 

TEMPRIS 

0 

1 

- Temp's from Individual Ply Heat Generation 

DAMP OUT 

0 

1 

- Suspends Damping's Output Printout 

FREQ 

STRAINS 

50.0 

0.0 

0.0 

0.0 

CURVTUR 

0.0 

0.0 

0.0 

1 $ Coefficient of 

friction 

is equal to 0.30 and the percentage of 

5 broken 

fibers is equal 

to 20 %. 

PLY 

1 

0.30 

0.20 

PLY 

2 

0.30 

0.20 

$ Axial stress amplitude of 20 ksi, (In-plane Loads Only). 

FMEMB 

400.0 

00. 

0 . 00.0 

PBEND 

0.0 

0.0 

o 

o 

o 

o 

PTRAN 

$ 

0.0 

0.0 

o 

o 

o 

PRINT 

ALL 

Output 

request for complete results (OPTION 0) 


Problem 2 

This problem includes the temperature rise option based on the individual strain 
energy dissipation per ply. The off-axis composite is a HT-S fiber/IMIIS matrix with a 
0.50 fiber volume ratio, subjected to a continuous sinusoidal excitation of a Qxx equal to 
10 ksi and a frequency of 50 Hz. The thickness of the unidirectional composite is 
0.04 in., hence, N x equals 400 lb/in., and the fiber orientation angle, 9^ equals 60 . 
Friction effects from broken fibers are ignored (FR1CTN = 0). The input data file is 
listed below: 
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s 

$ The temperature rise based on the individual strain energy 
$ dissipation per ply is considered in this example. The temperature 
$ difference is measured versus the room temperature of 70 degrees F. 


Demonstration Problem No. 2 
T 


I CAN /DAMP : 

COMSAT 

C3ANB F 

BIDE F 

RINDV F 

NONUDF T 

DEFECT F 

$ The composite is 0.04 inches thick and the fiber orientation 

$ angle is equal to 60 degrees. 

PLY 1 1 70.0 70.0 0.0 

p L Y 2 1 70.0 70.0 0.0 

$ HT-S fiber in Intermediate Modulus High Strength matrix ( IMHS ) . 

1HT-SIMHS 0.50 0.OOT3OOIMHS 0.00 0.00 

0 IF: 1 - Include Friction Effects 


60.0 

60.0 


0.02 

0.02 


MAT CRD 
FRIC'TN 
TEMPRIS 
DAMPOUT 
$ Frequency 


0.00 


Is equal to 50 Hz. 


1 - Temp ' s from Individual Ply Heat Generation 
1 - Suspends Damping's Output Printout 



FREQ 

STRAINS 

50.0 

0.0 

0.0 

0.0 


CURVTUR 

0.0 

0.0 

0.0 


PLY 

1 

0.00 

0.00 

. 

PLY 

2 

0.00 

0.00 


$ Axial 

stress amplitude 

of 10 ksi 

_ 

PMEMB 

400.0 

00. 

0.0 


PBEND 

0.0 

0.0 

0.0 


PTRAN 

0.0 

0.0 

o. 


S 

PRINT 

ALL 

Output 

request 


0.0 

0.0 

00 . 


for complete results (OPTION 0) 


Problem 3 


The variations of temperature and moisture will influence the damping performance of 
composite laminates. These hygrothermal effects on the flexural SDC of a 0.50 fiber- 
volume-ratio, HMSF/DX209 resin, quasi-isotropic (40/130/85/— 5) g laminate are modeled 
as follows: The temperature difference is selected as 200 F above reference temperature 
(70 F) by assigning a use temperature of 270 F while the cure temperature remains 
constant. The amount of moisture percent in the composite is 2 percent by weight. The 
laminate is subjected to an in-plane bending moment (free flexure) of M, = 2 lb m./in., 
and the thickness of each ply is 0.01 in. This example is illustrated in the following 

input data set: 


! 


| 

i 

i 
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$ 




$ This problem illustrates the hygrothermal effects on the flexural 

$ SDC of a 

(TH/90+TH/ 45+TH/-45+TH ) s laminate- The temperature differ- 

$ ence is 
$ 

equal to 200 F and the 

moistuie content is equal to 2 % . 

ICAN/DAMP: 

Demonstration Problem 

No. 3 

COMSAT 

T 



CSANB 

F 



BIDE 

F 



RINDV 

F 



WONUDF 

T 



DEFECT 

F 



$ Each ply 

is 0.01 

inches thick and the fiber orientation angle is 

$ equal to 

40 degrees. 


PLY 

1 

1 270.0 

70.0 2.0 40.0 0.01 

PLY 

2 

1 270.0 

70.0 2.0 130.0 0.01. 

PLY 

3 

1 270.0 

70.0 2.0 85.0 0.01 

PLY 

4 

1 270.0 

70.0 2.0 -5.0 0.01 

PLY 

5 

1 270.0 

70.0 2.0 -5.0 0.01 

PLY 

6 

1 270.0 

70.0 2.0 85.0 0.01 

PLY 

7 

1 270.0 

70.0 2.0 130.0 0.01 

PLY 

8 

1 270.0 

70.0 2.0 40.0 0.01 

1 $ HMSF fiber in DX209 epoxy resin 

(0.50 fiber volume ratio). 

MATCRD 

1HMSFD209 0.50 

0 . 00T300IMHS 0.00 0.00 0.00 

FRICTN 

0 

IF: 1 - Include Friction Effects 

TEMPRIS 

0 

1 - Temp r 

s from Individual Ply Heat Generation 

DAMPOUT 

0 

1 - Suspends Damping's Output Printout 

FREQ 

50.0 



STRAINS 

0.0 

0.0 0.0 


CURVTUR 

0.0 

0.0 0.0 


PLY 

1 

0.00 0.00 


PLY 

2 

0.00 0.00 


PLY 

3 

0.00 0.00 


PLY 

4 

0.00 0.00 


PLY 

5 

0.00 0.00 


PLY 

6 

0.00 0.00 


PLY 

7 

0.00 0.00 


PLY 

8 

0.00 0-00 


$ In-plane 

bending 

moment of 2 lb 

in/in. 

PMEMB 

0.0 

0.0 0.0 

0.0 

PBEND 

2.0 

0.0 0.0 

0.0 

FTRAN 

$ 

0.0 

0.0 0 . 

00 . 

PRINT 

ALL 

Output request 

for complete results (OPTION 0) 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, May 14, 1992 




Appendix A. — Symbols 


[A d UA] 

[C D ],[C] 

[D d ],[D] 

c 

E 

[E c ] 

G 

r 

K 

k 

Nr 


fb 


Ni 

S 

T 

t 

W,<5W 

z 

a 


4 

Subscripts 

b 

C 

c 

D 

f 


3x3 extensional damping and stiffness matrices, respectively 

3x3 coupling damping and stiffness matrices, respectively 

3x3 flexural damping and stiffness matrices, respectively 

thermal capacity 

normal modulus 

off-axis ply stiffness matrix 

shear modulus 

forcing frequency, Hz 

thermal conductivity 

volume ratio 

fraction of broken fibers in a ply 

number of plies 

strength 

temperature 

ply thickness 

stored and dissipated strain energies 

distance from laminate midplane surface 

thermal expansion coefficient 

engineering strain 

midplane engineering strain vector 

curvature vector 

fiber- matrix friction coefficient 

Poisson’s ratio 

mass density 

stress 

specific damping capacity 

ply bottom 

compression 

(ply) off-axis composite 

damping 

fibers 


1G 


Ill 

n 

o 

s 

t 

x,y,z 
xx, yy 
SS 

1 , 2,3 

11 

22 

33 

12 

23,13 


(ply) on-axis composite 
matrix 

normal direction 
reference 
shear direction 
ply top 

structural axes 

in-plane normal, off-axis, x and y directions, respectively 
in-plane shear, off-axis 
material axes 

normal longitudinal direction 
normal in-plane transverse direction 
normal out-of-plane transverse direction 
shear in-plane direction 
shear out-of-plane directions 
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Appendix B. — Sample ICAN/DAMP Input Data Set 


$ Title Card Card Group I * 

S 

ICAN/DAMP: QUASI -ISOTROPIC LAMINATE , 8 ' PLIES , DT=000 F,M=0% 

$ 

5 Boolean Card Group Card Group II* 

$ 

COMSAT T 

CSANB F 

$ The logical variable BIDE is Not being used, {Interply Layer 
$ Distortion Energy Coefficient (ILDC = 0)}. 

BIDE F 

RINDV F 

NONUDF T 

DEFECT F 

$ 

5 Laminate Configuration Card Group Card Group III. 

$ The following is the PLY card group. The laminate configuration is 
S specified as follows: 


70.0 

0.0 

00.0 

0.01 

70.0 

0.0 

90.0 

0.01 

70.0 

0.0 

45.0 

0.01 

70.0 

0.0 

-45.0 

0.01 

70.0 

0.0 

-45.0 

0.01 

70.0 

0.0 

45.0 

0.01 

70.0 

0.0 

90.0 

0.01 

70.0 

0.0 

00.0 

0.01 

■ Card 
)e used 

Group IV 
in the 

analys is 

are 


described as follows: 

MATCFD 1HMSFD210 0.50 0.00T300IMHS 


0.00 0.00 0.00 


5 

£ ************************************* A************** *************** 

$ *** Damping Analysis Card Groups. *** 

$ The Input data cards needed to perform the damping analysis are 
$ specified below; 

$ 

5 Damping Option Card Group Card Group V. 


FRICTN 

TEMPRIS 

DAKPOUT 


IF: 1 - Include Friction Effects 

1 - Temp ' s from Individual Ply Heat Generation 
1 - Suspende Damping's Output Printout 


$ Forcing Frequency Card Group VI. 

$ 

FREQ 50.0 

$ 

$ Deformation Card Group VII. 


STRAINS 

0.0 

0.0 

0.0 

CURVTUR 

0.0 

0.0 

0.0 

$ 

$ Friction Data 

r {MU 

and NFE 

$ 

PLY 

1 

0.30 

0.20 

PLY 

2 

0.30 

0.20 

PLY 

3 

0.30 

0.20 

PLY 

4 

0.30 

0.20 

PLY 

5 

0.30 

0.20 

PLY 

6 

0.30 

0.20 

PLY 

7 

0.30 

0.20 

PLY 

8 

0.30 

0.20 


* 

$ *** END of Damping Analysis Card Groups. *** 

^ ********************************* *********************** *********** 

$ 

$ Loads Card Group Card Group IX. 
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S In-plane Membrane Loads. 

PMEMB 00.0 00. 0.0 0.0 

$ In-plane Bending Loads. 

PBEND 2.0 0.0 0.0 0.0 

$ Transverse Loads. 

PTRAN 0.0 0.0 0. 00. 

$ 

5 Cyclic or Fatigue Load Data Group Card Group X. 

$ Loading for the Durability Analysis. Meehan ical /Thermal Fatigue. 


CNXX 

200. 

100. 

100. 

0.1 

CNYY 

-50. 

-100. 

10. 

0.1 

C'NXY 

20. 

10. 

100. 

0.2 

CMXX 

10. 

5. 

10. 

0.01 

CMYY 

1 . 

2. 

1000. 

0.15 

CMXY 

2. 

1 . 

100. 

0.01 


$ 

5 Output Selection Data Group Card Group XI. 

$ The output can be tailor. OPTION 0 is the output request for 

$ complete results. 

5 

PRINT ALL Output request for complete results {OPTION 0) 
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Appendix C. — Sample Data Bank 


i 

j 


FIBER PROPERTIES 
T300 GRAPHITE FIBER. 

$ 

$ 

$ 

Humber of fibers per end 
Filament equivalent diameter 
Weight density 
Normal moduli {11) 

Normal moduli {22) 

Poisson's ratio (12) 
Poisson's ratio (23) 

Shear moduli (12) 

Shear moduli (23) 

Thermal expansion coef. (11) 
Thermal expansion coef. (22) 
Heat conductivity (11) 

Heat conductivity (22) 

Heat capacity 
Fiber tensile strength 
Fiber compressive strength 
Fiber damping capacity (11) 
Fiber damping capacity (22) 
Fiber damping capacity (12) 
Fiber damping capacity (23) 


Nf 

3000 

df 

0. 300E-03 

Rliof 

0 . 640E-01 

Ef 11 

0. 320E+08 

Ef 22 

0 . 200E + 07 

N u f 1 2 

0.200E+00 

Muf 23 

0. 250E4-00 

Gf 12 

0. 130E+07 

Gf23 

0 . 700E ^6 

Alfafll 

- . 550E-06 

Alfaf22 

0 . 560E-05 

Kf 11 

0. 403E + 01 

Kf 22 

0.403E+00 

Cf 

0.170E+00 

SfT 

0. 350E^06 

SfC 

0. 300E+06 

Zetaf 11 

0. 100E-02 

Zetaf 22 

0. 000E+00 

Zetaf 12 

0 . 000E+00 

Zetaf23 

O.OOOE+OO 


number 
inches 
lb/ in* * 3 
psi 
psi 

non-dim 
non-dim 
ps i 
psi 

in/in/F 

in/in/F 

BTU-in/hr/J n**2/F 

BTU-in/hr/ in* *2/F 

BTU/lb/F 

psi 

psi 

non-dim 

non-dim 

non-dim 

non-dim 


I 

i 

j 


HMSF HIGH MODULUS SURFACE TREATED FIBER. 

S 

5 


S 

Number of fibers per end 
Filament equivalent diameter 
Weight density 
Normal moduli (11) 

Normal moduli (22) 

Poisson's ratio (12) 
Poisson's ratio (23) 

Shear moduli (12) 

Shear moduli (23) 

Thermal expansion coef. (11) 
Thermal expansion coef. (22) 
Heat conductivity (11) 

Heat conductivity (22) 

Heat capacity 
Fiber tensile strength 
Fiber compressive strength 
Fiber damping capacity (11) 
Fiber damping capacity (22) 
Fiber damping capacity (12) 
Fiber damping capacity (23) 


Nf 

10000 

number 


df 

0 . 300E-03 

inches 


Rhof 

0 . 703E-01 

lb/iu* * 3 


Ef 11 

0 . 550E+08 

psi 


E f 22 

0 . 900E+06 

ps L 

_ 

Nuf 12 

0. 200E+00 

non -dim 


Nuf23 

0.250E+00 

non-dim 

* 

Gf 12 

0. 110E + 07 

ps i 

« 

Gf 2 3 

0.700E+06 

psi 

- 

Alfafll 

- . 550E-06 

in/in/F 


Alfaf22 

0 . 560E-05 

in/in/F 

“ 

Kf 11 

0 . 403E *01 

BTU-in/hr/in* *2/F . 


Kf 22 

0. 403E+00 

BTU-in/hr/in** 2 /F 

= 

" 

Cf 

0. 170E+00 

BTU/lb/F 


SfT 

0 . 280E+06 

psi 


SfC 

0.200E+06 

psi 


Zetaf 11 

0 . 400E-02 

non-dira 

= 

Zetaf22 

0.400E-02 

non-dim 


Zetaf 12 

0. 400E-02 

non-dim 


Zetaf 23 

0. 400E-02 

non-dim 



HT-S 

$ 

S 

$ 

Number of fibers per end 
Filament equivalent diameter 
Weight density 
Normal moduli (11) 

Normal moduli (22) 

Poisson's ratio (12) 
Poisson's ratio (23) 

Shear moduli (12) 

Shear moduli (23) 

Thermal expansion coef, (11) 
Thermal expansion coef. (22) 
Heat conductivity (11) 


Nf 

10000 

number 


df 

0. 300E-03 

inches 


Rhof 

0 . 640E-01 

lb/ in* * 3 


Ef 11 

0. 320E + 08 

ps i 


Ef 22 

0 , 250E + 07 

ps i 

_ 

Nuf 12 

0 . 200E + 00 

non-dim 


Nuf 23 

0 . 250E+ 00 

non-dim 

— 

Gf 12 

0. 150E+07 

psi 


Gf 23 

0. 100E+07 

psi 


Alfafll 

-. 550E-06 

in/in/F 


Al f af 22 

0 . 560E-05 

in f in/F 

- 

Kf 1 1 

0 . 403E+01 

BTU-in/hr /in* * 2 /F 
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Heat 

conductivity (22) 


Kf 22 

0 . 403E+00 

BTU-in/hr/ in* * 2/F 

Heat 

capacl ty 



Cf 

0 . 17 0E + 00 

BTU/lb/F 

Fiber 

tensile 

strength 


SfT 

0 . 350E+ 06 

psi 

Fiber- 

compressive strength 

SfC 

0 . 300E+06 

psi 

Fiber 

damping 

capacity 

<H) 

Zetafll 

0.700E-02 

non-dim 

Fiber 

damping 

capacity 

(22) 

Zetaf22 

0 . 700E-02 

non -dim 

Fiber 

damping 

capacity 

(12) 

Zetafl2 

0. 600E-02 

non-dim 

Fiber 

damping 

capacity 

(23) 

Zetaf23 

0. 600E-02 

non-dim 

KEVL 

KEVLAR 

FIBER. 






$ 

$ 

s 


Number of fibers per end 

Nf 

580 

number 

Filament equivalent diameter 

df 

0 . 460E-03 

Inches 

Weight density 


Rhof 

0 . 530E-01 

lb /in* *3 

Normal moduli (11) 


Ef 11 

0. 220E* 08 

psi 

Normal moduli (22) 


Ef 22 

0 . 600E+06 

psi 

Poisson's ratio (12) 


Nuf 12 

0.350E+00 

non -dim 

Poisson's ratio (23) 


Nuf 23 

0 . 350E+00 

non-dim 

Shear moduli (12) 


Gf 12 

0.420E+06 

psi 

Shear moduli (23) 


Gf 23 

0 . 220E * 06 

psi 

Thermal expansion coef. 

(11) 

Alfafll 

- . 220E-05 

in/in/F 

Thermal expansion coef. 

(22) 

A1 faf 22 

0 . 300E-04 

in '’in/F 

Heat conductivity (11) 


Kf 11 

0. 118E-01 

BTU-in/hr/in* *2/F 

Heat conductivity (22) 


Kf 22 

0. 118E-01 

BTU-in/hr /in* *2/F 

Heat capacity 


Cf 

0.250E+00 

BTU/lb/F 

Fiber tensile strength 


SfT 

0. 400Et06 

ps i 

Fiber compressive strength 

SfC 

0.750Ef05 

ps i 

Fiber damping capacity 

(ID 

Zetafll 

0. 180E-01 

non-dim 

Fiber damping capacity 

(22) 

Zetaf 22 

0. 180E-01 

non-dim 

Fiber damping capacity 

(12) 

Zetaf 12 

0. 180E-01 

non-dim 

Fiber damping capacity 

(23) 

Zetaf23 

0- 180E-01 

non-dim 


EGLA E-GLASS FIBER. 
$ 

$ 

$ 


Number of fibers per end 

Nf 

204 

number 

Filament equivalent diameter 

df 

0 . 360E-03 

inches 

Weight density 


Rhof 

0 . 900E-01 

1 b / i n * * 3 

Normal moduli (11) 


Ef 1 1 

0 . 105E ♦ 00 

ps i 

Normal moduli (22) 


Ef 22 

0. 105E+08 

psi 

Poisson's ratio (12) 


Nuf 12 

0 . 200E+00 

non-dim 

Poisson's ratio (23) 


Nuf 2 3 

0 . 200E+00 

non-dim 

Shear moduli (12) 


Gf 12 

0.437E+07 

psi 

Shear moduli (23) 


Gf23 

0 . 437E + 07 

psi 

Thermal expansion coef. 

(11) 

Alfafll 

0 . 280E-05 

in/in/F 

Thermal expansion coef. 

(22) 

Alfa f 22 

0 - 280E-05 

in/in/F 

Heat conductivity (11) 


Kf 11 

0.521E-01 

BTU-in/hr/ in* * 2/F 

Heat conductivity (22) 


Kf 22 

0.5 2 IE-01 

BTU-in/hr/ in* * 2 /F 

Heat capacJ ty 


Cf 

0. 170E+00 

BTU/lb/F 

Fiber tensile strength 


SfT 

0. 360E+0G 

psi 

Fiber compressive strength 

SfC 

0.360E406 

psi 

Fiber damping capacity 

(11) 

Zetafll 

0. 110E-01 

non -dim 

Fiber damping capacity 

(22) 

Zetaf22 

0 . 11QE-01 

non-dim 

Fiber damping capacity 

(12) 

Zetaf 12 

0. 600E-02 

non-dim 

Fiber damping capacity 

(23) 

Zetaf 23 

0. 600E-02 

non-dim 


OVER END OF FIBER PROPERTIES 


MATRIX PROPERTIES 

LMLS LOW MODULUS LOW STRENGTH MATRIX. 
$ 

$ 

$ 

Weight density Rhom 

Normal modulus Em 

Poisson's ratio Num 


0 . 420E-01 
0 . 320E+06 
0 . 4 30E+-00 


lb/in* * 3 
psi 

non-dim 



Thermal expansion coef. 
Matrix heat conductivity 
Heat capacity 
Matrix tensile strength 
Matrix compressive strength 
Matrix shear strength 
Allowable tensile strain 
Allowable compr. strain 
Allowable shear strain 
Allowable torsional strain 
Void heat conductivity 
Glass transition temperature 
Matrix damp. Capacity-Formal 
Matrix damp. capacity-Shear 


Alfa ra 

0 . 570E-04 

Km 

0 . 868E-02 

Cm 

0. 250E+00 

SmT 

0 . 800E+04 

SmC 

0. 150E+05 

SmS 

0 . 800E+04 

eps mT 

0 . 810E-01 

eps mC 

0. 150E + Q0 

eps mS 

0. 100E+00 

eps mTOR 

0. 100E+00 

kv 

0.225E+00 

Tgdr 

0. 350E+03 

Zeta mn 

0. 850E-01 

Zeta ms 

0.138E+00 


in/in/F 

BTU-in/hr/in* * 2/F 

BTU/lb/F 

psi 

psi 

psi 

in/in 

in /in 

in/in 

in /in 

BTU- in/hr/ in* * 2/F 
F 

non-dim 

non-dim 


I 

! 


IMLS 

$ 

$ 


INTER. MODULUS LOW STRENGTH (SIMILAR TO 350 B) 


•? 

Weight density 

Rhom 

0. 460E-01 

lb/in* *3 

Normal modulus 

Em 

0 . 500E + 06 

psi 

Poisson's ratio 

Num 

0. 410E+00 

non-dim 

Thermal expansion coef. 

Alfa in 

0.570E-04 

in/in/F 

Matrix heat conductivity 

Km 

0. 868E-02 

BTU-in/hr/in* *2/F 

Heat capacity 

Cm 

0.250E+00 

BTU/lb/F 

Matrix tensile strength 

SmT 

0.700E+04 

psi 

Matrix compressive strength 

SmC 

0.210E+05 

ps L 

Matrix shear strength 

SmS 

0 . 700E+ 04 

psi 

Allowable tensile strain 

eps mT 

0. 140E-01 

in/in 

Allowable compr, strain 

eps mC 

0 . 420E-01 

in/in 

Allowable shear strain 

eps mS 

0 . 320E-01 

in/in 

Allowable torsional strain 

eps mTOR 

0 . 320E-01 

in/in 

Void heat conductivity 

kv 

0.225E+00 

BTU-in/hr/in* * 2/F 

Glass transition temperature 

Tgdr 

0 . 420E+03 

F 

Matrix damp, capacity -Formal 

Zeta mn 

0. 850E-01 

non-dim 

Matrix damp. capacity-Shear 

Zeta ms 

0. 138E+00 

non-dim 

XMHS INTERMEDIATE MODULUS HIGH STRENGTH 
$ 

MATRIX. 


5 

Weight density 

Rhom 

0. 440E-01 

Ib/in* *3 

Normal modulus 

Em 

0 . 500E+06 

psi 

Poisson's ratio 

Num 

0. 350E+00 

non-dim 

Thermal expansion coef. 

Alfa m 

0. 360E-04 

in/ in/F 

Matrix heat conductivity 

Km 

0 . 868E-02 

BTU- in/hr/ in* * 2/F 

Heat capacity 

Cm 

0 . 250E+00 

BTU/lb/F 

Matrix tensile strength 

SmT 

0. 150E + 05 

psi 

Matrix compressive strength 

SmC 

0. 350E+05 

psi 

Matrix shear strength 

SmS 

0. 130E+05 

psi 

Allowable tensile strain 

eps mT 

0 . 200E-01 

in/in 

Allowable compr. strain 

eps mC 

0 . 500E-01 

in/in 

Allowable shear strain 

eps mS 

0. 350E-01 

in/in 

Allowable torsional strain 

eps mTOR 

0. 350E-01 

in /in 

Void heat conductivity 

kv 

0 . 225E+00 

BTU-in/hr/in* * 2/F 

Glass transition temperature 

Tgdr 

0. 420E+03 

F 

Matrix damp. capacity-Formal 

Zeta mn 

0 . 660E-01 

non-dim 

Matrix damp, capacity -3 hear 

Zeta ms 

0.690E-01 

non-dim 


K 

4 


D209 DX209 RESIN 

$ 




$ 

S 

Weight density 

Rhom 

0. 440E-01 

lb/in* *3 

Normal modulus 

Em 

0 . 500E+06 

psi 

Poisson's ratio 

Num 

0 . 350E + 00 

non- dim 

Thermal expansion coef. 

Alfa m 

0 . 360E-04 

in/in/F 

Matrix heat conductivity 

Km 

0 . 868E-02 

BTU-in/hr/in** 2/F 

Heat capacity 

Cm 

0 . 250E+00 

BTU/lb/F 

Matrix tensile strength 

SmT 

0. 150E+05 

psi 
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Matrix compressive strength 

SmC 

0. 350E+05 

psi 

Matrix shear strength 

SmS 

0. 1 30E+05 

psi 

Allowable tensile strain 

eps mT 

0 . 200E-0 1 

in/in 

Allowable compr. strain 

eps mC 

0 . 500E-01 

in/in 

Allowable shear strain 

eps mS 

0 . 350E-01 

in/in 

Allowable torsional strain 

eps mTOR 

0. 350E-01 

in /in 

Void heat conductivity 

kv 

0. 225E+00 

BTU- in/hr /in * * 2 /F 

Glass transition temperature 

Tgdr 

0.420E+03 

F 

Matrix damp. capacity-Formal 

Zeta mn 

0.172E+00 

non-dim 

Matrix damp, capacity-Shear 

Zeta ms 

0.168E+00 

non-dim 


D210 DX210 RESIN 

5 
$ 

$ 


Weight density 

Rhom 

0. 440E-01 

lb/in* * 3 

Normal modulus 

Em 

0.500E+06 

psi 

Poisson's ratio 

Num 

0 . 350E + 00 

non-dim 

Thermal expansion coef. 

Alfa m 

0 . 360E-04 

in/in/F 

Matrix heat conductivity 

Km 

0 . 868E-02 

BTU- In /hr /In* * 2/F 

Heat capacity 

Cm 

0. 250E+00 

BTU/lb/F 

Matrix tensile strength 

SirtT 

0. 150E+05 

ps i 

Matrix compressive strength 

SmC 

0. 350E+05 

psi 

Matrix shear strength 

Sm3 

0. 130E+05 

psi 

Allowable tensile strain 

eps mT 

0. 200E-01 

in/in 

Allowable compr. strain 

eps mC 

0 . 500E-0 1 

in/in 

Allowable shear strain 

eps mS 

0 . 350E-01 

in /in 

Allowable torsional strain 

eps mTOR 

0 . 350E-0 1 

in/in 

Void heat conductivity 

kv 

0.225E+00 

BTU-in/hr/in* *2/F 

Glass transition temperature 

Tgdr 

0. 420E+03 

F 

Matrix damp. capacity-Formal 

Zeta mn 

0. 103E+00 

non-dim 

Matrix damp, capacity-shear 

Zeta ms 

0. 1I8E+00 

non -dim 


OVER END OF MATRIX PROPERTIES 
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Appendix D. — Sample Damping Module Output 


**************^*****M******* 

DAMPING A N A L Y 

A**************************** 


h******hmm»m***********m******h*************** 

SIS OF COMPOSITE LAMINATES 
***************************** ************************ 


SUPRESSED OUTPUT OPTIONS: 


- FRICTION EFFECT ON ZETAL11 

- TEMPERATURE PROFILES BASED ON INDIVIDUAL PLY DISSIPATION 


** ON-AXIS DAMPING ** 


i 

I 


MATERIAL SYSTEM 1: (HMSF/D210 @ DT= 0.00, M-0.00) 


FIBER DAMPING CAPACITY ( HMSF ) : 


ZFU = 

0. 4000E-02 

ZF12 = 

0 . 4000E-02 

ZF13 = 

0. 4000E-02 

ZF22 - 

0. 4000E-02 

ZF23 - 

0. 4000E-02 

ZF33 = 

0. 4000E-02 


MATRIX DAMPING CAPACITY (D210) f 


ZMN = 0 . 1030E + 00 

ZMS = 0 . 1 180E-* 00 


i 

I 

i 

i 

i 

i 


i 


| 

: COMPOSITE DAMPING CAPACITY (W/O FRICTION): 


a 


3 

ZL11 

= 

0 . 4892E-02 

S 

ZL12 

= 

0. 704 IE-01 


ZL13 

= 

0.7041E-01 


ZL21 

* 

0. 704 IE-01 

- 

ZL22 

- 

0 . 4202E-01 


ZL23 

= 

0 . 4817E-01 


ZL31 


0. 7041E-01 


ZL32 

= 

0. 4817E-01 

■ 

ZL33 

= 

0 . 4202E-01 
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ini in ii : i i mu ii i II I 1 1 III Mill Hit 



OFF-AXIS 


DAMPING ** 


* * 



PLY If 1 

PLY # 2 

FLY # 3 

PLY # 4 


HMSF/D2 10 

HMSF/D210 

HMSF/D2 10 

HMSF/D210 


0.0 DEG 

90.0 DEG 

45.0 DEG 

-45.0 DEG 

DAMPING MATRIX: ZCXX 

0.4892E-02 

0. 4202E-01 

0. 4693E-01 

0 . 4693E-01 

ZCXY 

0 . OOOOE+OO 

1509E-12 

2347E-01 

- . 2 347E-01 

ZCXS 

O.OOOOEkOO 

0 . 3598E-07 

- . 9283E-02 

0 . 9283E-02 

ZCYX 

0 . OOOOE+OO 

- . 1509E-12 

- . 2 347E-01 

2347E-01 

ZCYY 

0 . 4202E-01 

0. 4892E-02 

0 . 4693E-01 

0. 4693E-01 

ZCYS 

0 . OOOOE+OO 

8304E-07 

9283E-02 

0. 9283E-02 

zcsx 

0 . 0000E+ 00 

0 . 7 195E-07 

- . 1857E-01 

0. 1857E-01 

ZC'SY 

0 . 0000E+00 

- . 1661E-06 

1857E-01 

0. 1857E-01 

zcss 

0. 704 IE-01 

0 . 7041E-03 

0. 2346E-01 

0 . 2346E-01 

AXIAL DAMPING: ZCAXLX 

0 . 4892E-02 

0. 4202E-01 

0. 5923E-01 

0. 5923E-01 

ZC AXL 1 

0. 4202E-01 

0 . 4892E-02 

0. 5923E-01 

0. 5923E-01 

ZCAXLS 

0. 7041E-01 

0 . 704 IE-01 

0. 4083E-01 

0 . 4083E-01 



PLY # 5 

PLY # 6 

FLY If 7 

PLY # 0 


HMSF/D2 10 

HMSF/D210 

IIMSF/D2 10 

HMSF/D210 


-45.0 DEG 

45.0 DEG 

90.0 DEG 

0.0 DEG 

DAMPING MATRIX: ZCXX 

0 . 4693E-01 

0. 4693E-01 

0. 4202E-01 

0 . 4892E-02 

ZCXY 

2347E-01 

-. 2347E-01 

-. 1509E-12 

0. OOOOE+OO 

ZCXS 

0 . 9283E-02 

- . 9283E-02 

0. 3593E-07 

0 . OOOOE+OO 

ZCYX 

-.2347E-01 

- . 2347E-01 

-. 1509E-12 

0. OOOOE+OO 

ZCYY 

0. 4693E-01 

0. 4693E-01 

0 . 4892E-02 

0 . 4202E-01 

ZCYS 

0 . 9283E-02 

- . 9283E-02 

- . 03O4E-O7 

0. OOOOE+OO 

ZCSX 

0 . 1057E-O1 

1857E-01 

0 . 7195E-07 

0. OOOOE+OO 

ZCSY 

0.1857E-01 

1857E-01 

1661E-06 

0. OOOOE+OO 

ZCSS 

0 . 2346E-01 

0 . 2346E-01 

0.704LE-01 

0 . 704 IE-01 

AXIAL DAMPING: ZCAXLX 

0 . 5923E-01 

0 . 5923E-01 

0 . 4202E-01 

0 . 4892E-02 

ZCAXLY 

0 . 5923E-01 

0 . 5923E-01 

0 . 4092E-O2 

0. 4202E-01 

ZCAXLS 

0. 4063E-01 

0. 4083E-01 

0 . 7041E-01 

0. 7041E-01 
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LAMINATE DAMPING 


* * 


> LAMINATE MECHANICAL STRAINS < 


MIDPLANE 

STRAIN , {EO} 

CURVATURE, {K} 

E01 * 

0. 1989E-10 

El 

= 0.2777E-02 

E02 = 

-0.2072E-10 

K2 

= -0 . 2 107E-03 

E03 = 

0. 1596E-10 

K3 

= -0 . 1279E-02 


SPECIFIC DAMPING CAPACITY = 0. 60160160 IE-02 
FOR ABOVE STRAINS. 

DISSIPATED STRAIN ENERGY = 0 . 16704 1 14 4E-04 
(PER AREA-CYCLE) 

STORED STRAIN ENERGY = 0. 277660601E-02 

(PER AREA-CYCLE) 


> LAMINATE DAMPING MATRICES < 


EXTENSIONAL DAMPING: {ACD} COUPLING DAMPING: {BCD} 


0 . 6362E+04 
0 . 6840E+ 03 
0. 1009E-02 

0 . 6840E+03 
0 . 6362E+04 
0. 1642E-02 

0 . 8 377E-03 
0. 1623E-02 
0. 2839E+04 

0. 1526E-04 
0.441 IE-05 
0 . 6083E-06 

0 . 4768E-05 
0 . 1526E-04 
0 . 7982E-06 

0. 1514E-06 
0. 3414E-06 
0. 4005E-04 

COUPLING DAMPING: 

{CCD} 

FLEXURAL DAMPING: 

{DCD } 

0. 1526E-04 
0.441 IE-05 
0 . 6083E-06 

0. 4768E-05 
0. 1526E-04 
0 . 7982E-06 

0. 1514E-06 
0. 34 14E-06 
0. 4005E-04 

0 . 4 148E+ 01 
0. 1973E + 00 
0. 1128E + 00 

0 . 2868E+00 
0 . 2884E + 01 
0. 1128EfOO 

0 . 9786E-01 
0 . 9786E-01 
0. 1391E+01 
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** PLY 

T E M P E 

RATURE DIS 

T R I 8 U T 

ION ** 

FREQUENCY 

(HZ) = 

50.00 



PLY 

INITIAL 

T - TO 

T - TO 

T - TO 


TEMP (TO) 

(AVERAGE) 

(BOTTOM) 

(TOP) 


(DEG. F) 

(DEG. F) 

( DEG . F ) 

(DEG. F) 

1 

70.000 

0.000 

0.000 

0.000 

2 

70.000 

0.001 

0.000 

0.001 

3 

70.000 

0.001 

0.001 

0.001 

4 

70.000 

0.001 

0.001 

0.001 

5 

70.000 

0.001 

0.001 

0.001 

6 

70.000 

0.001 

0.001 

0.001 

7 

70.000 

0.001 

0.001 

0.000 

8 

70.000 

0.000 

0.000 

0.000 

:***HH*********************A*****M**HH***************M 

************** 

END 0 

F DAM 

PING ANALY 

SIS 


:**HHH****** 

*********** 

*********************, 

************* 

************** 
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Appendix E. — Files for Compiling and Executing 
ICAN/DAMP Under VM/SP Operating System 
and VAX /VMS Environment 

*REXX EXEC * file for compiling ICAN/DAMP under 
VM/SP Operating System. 


/* COMPILE EXEC */ 

ARG pliylen phytyp phytuod junk 
'set eras type lit 

F0RTVS2 phylen '{MAP FLAG (E) opt (2) NOPRINT' 
rcode - re 

if rcode = 4 then rcode = 0 
'set crastype it' 
load phylen '(clear' 
genraod phylen 
erase phylen text phymod 
exit rcode 


'RSXX EXEC* file for execut i ng ICAN/DAMP under 
VM/SP Ope ra ting System. 


1 








*/ 


/* 

x — 








/* 

X 





X 

*/ 

: 

/* 

X 



ICAN/DAMP EXEC 

X 

*/ 

i 

/* 

X 





X 

*/ 

55 

/* 

X 

This 

Is a 

REXX 

EXEC file to run ICAN/DAMP 

X 

*/ 


/* 

X 

File 

unit 

7 is 

for INHYD - 

X 

*/ 

- 

/* 

X 

File 

unit 

5 is 

for Input. 

X 

*/ 

= 

/* 

X 

File 

unit 

8 is 

for Data Bank. 

X 

*/ 


/* 

X 

File 

unit 

6 is 

for Output. 

X 

*/ 


/* 

X 




filename: 'iendamp'. 

X 

*/ 

- 

/* 

X 

Fortran 

code 

X 

*/ 

- 

/* 

X 


Data 

bank 

filename: 'tbledam*. 

X 

*/ 


/* 

X 





X 

*/ 


/* 

X 

NOTE: 

: The 

data 

bank filename Is specified 

X 

*/ 


/* 

X 


In subroutine OPENFL. 

X 

*/ 


/* 

X— 






*/ 


parse ary fa fra ft fout . 

if fout = ' ' then fout = 'dampout ' 

'FILEDEF FT07F001 DISK JUNK DAT D1 (LRECL 132 XTENT 9999' 
'FILEDEF FT05F001 DISK 'fn' DATA Di' 

'FILEDEF FT06F001 DISK 'fout' LISTING Dl (LRECL 132 XTENT 9999' 
'iendamp ' 

'fil * clear' 

'ERASE JUNK DAT Dl' 
exit rc 


III 
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DCL Command Procedures File for executing ICAN/DAMP under VAX/VMS Operating System. 


$J — 

5! COM FILE FOR THE EXECUTION OF ICAN/DAMP ON VAX 

$1 

$J 

$1 Input files 

$ assign icninp.dat for005 

$ assign junk.dat for007 

$! NOTE: The data bank file is specified 

$! in subroutine OPENFL. 

$! Output files 

$ assign dampout.dat for006 

$ assign sys$output forOOl 

SJ run Icndmp.exe 

$ run icudinp 
$! 

$ delete junk.dat;* 

$ exit 
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Appendix F. — On- Axis Damping, Off-Axis Damping, and 
Laminate 

Damping Equations 


Damping Micromechanics 

Longitudinal normal ply SDC 'I’m'- 


</>! 


II 



E fll 

E U 1 


E„ 


+ tp 


L 

inn m 


E 


111 


( 1 ) 


Transverse normal ply SDC’s tp n2 , 0 133 : 


V’122 - ^(22 



^inn 



( 2 ) 


^133 = ^122 


( 3 ) 


where 


E 2 2 = (1 ) E m + 






1 - [1 - (E m /E f22 )] 


( 4 ) 


Tn-plane shear ply SDC 

I ^12 / ^12 

^112 " %2 V k f p + Ks ( 1 ~ V k f ) p— 


3 fI2 


( 5 ) 


where 


'12 


(1 ~ff~) 


G m + 




1 [1 -(G m /G fl2 )] 


( 6 ) 
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Out-of-plane shear ply SDC’s t/i| 23 , ^113 


I G 123 

^123 " ^f23 V k f 7 T- + * ms ( 1 
G f23 



( 7 ) 


^113 = ^112 


where 


G 123 


E 122 

2 ( 1 + "123 ) 


( 8 ) 


( 9 ) 


E, 


122 



(10) 


*'123 = 


1 ~ k f "in 


"f23 " 


( 1 - k f ) "m 
1 " k f"m 


( 11 ) 


In-Plane Damping 

Off-axis composite damping [0j: 

where 

l R 


W c \ = ( R <r ] T [^i)[R <r " 1 ) T 


cos^ 9 

sin^0 

sin 29 

sin^fl 

cos^ 9 

-sin 29 

-0.5 sin 2 0 

0.5 sin 2(9 

cos 2 9 


( 12 ) 


(13) 


On-axis damping matrix [ t/’i] 
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0111 

0 

0 

w = 

0 

0122 

0 


0 

0 

0(12 

hence, the off-axis damping matrix is: 


V’cxx 

{!) 

v cxy 

0C xs 

w = 

V’cyx 

V’cyy 

0cys 


^csx 

^csy 

0C9S 


( 14 ) 


( 15 ) 


Laminate Damping 

Extensional Damping Matrix [A D |: 

|a d 1 


Ni 

E K -*bhI E cliWi 

i=i 


Coupling Damping Matrix (C D ): 


|C D 1 


Nl 

i=i 2 


Flexural Damping Matrix |Dp|: 


! d d1 


N, 

E ;K 3 -«b 

i=l 3 


where 


( 16 ) 


( 17 ) 


( 18 ) 
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[E c ] 


-1 


, 1/ V 

1 cyx _ csx 

F F O 

^cxx ^cyy ^cxy 

v \ v 

cxy 1 csy 

F F P 

cxx cyy cxy 

v v , 

cxs cys 1 

F r n 

°cxx ^cyy ^cxy 


Dissipated Strain Energy (Per area-cycle) 6 W: 


1 rt, 

[a d ] 

|C D l' 


U 

SW = 2 e° K 





2 

jC D l 

I°dI 


[*) 


Stored Strain Energy (Per area-cycle) W: 


W = A e° k 

|A| 

|CJ 


■° 

2 

JC) 

ID], 


[«) 


Laminate SDC tp: 



(19) 


( 20 ) 


( 21 ) 


( 22 ) 
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Appendix G. — Damping Module Structure 


Subroutine DAMP is the main subroutine of the damping module and is called from 
subroutine ICANSB of the main 1C AN/D AMP code. A brief description of the module’s 
basic subroutines is given as follows: 

(1) Subroutine DBNKR2 — This subroutine reads the data base of constituent 
materials and retrieves the Tiber and matrix specific damping capacities. 

(2) Subroutine ORTDMP— Calculates the on-axis ply damping capacities from the 
basic composite material constituent properties and fiber content (micromechanics). 

(3) Subroutine LAMDMP— Calculates the off-axis SDC’s of each ply, the laminate 
damping matrices [A D ], [C D J, and [D d ], and the temperature rise (optional) based on the 
individual strain energy dissipation in each ply. 

(4) Subroutine LAMLOS — Calculates the laminate SDC for a given local deforma- 
tion and the temperature rise through-the-thickness based on the average strain energy 
loss. 


v 

v 
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